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ABSTRACT 

Planck's all-sky surveys at 30-857 GHz provide an unprecedented opportunity to follow the radio spectra of a large sample of extragalactic 
sources to frequencies 2-20 times higher than allowed by past, large-area, ground-based surveys. We combine the results of the Planck Early 
Release Compact Source Catalog (ERCSC) with quasi- simultaneous ground-based observations as well as archival data at frequencies below 
or overlapping Planck frequency bands, to validate the astrometry and photometry of the ERCSC radio sources and study the spectral features 
shown in this new frequency window opened by Planck. The ERCSC source positions and flux density scales are found to be consistent with the 
ground-based observations. We present and discuss the spectral energy distributions (SEDs) of a sample of "extreme" radio sources, to illustrate 
the richness of the ERCSC for the study of extragalactic radio sources. Variability is found to play a role in the unusual spectral features of some 
of these sources. 

Key words. Surveys: radio sources — Radio continuum: galaxies — Radiation mechanisms: general 

1, Introduction Results" Planck papers there are three that address the extra- 

galactic radio source population. Planck Collaboration (201 Id) 
This paper is one of a series based on observations of com- examines statistical properties such as number counts and spec- 
pact sources by the Planck satellite that are included in the ^^.^1 index distributions, but only at frequencies > 30 GHz. Planck 
ERCSC (Planck Collaboration 2011c). Among these "Early Collaboration (201 le) incorporates Planck measurements and 

supporting ground-based and satellite observations to refine 

* Corresponding author: B. Partridge, bpartrid@haverford.edu models for the physical properties of a sample of -100 bright 

^ Planck (http://www.esa.int/Planck) is a project of the blazars. Here we address the observed diverse, sometimes pe- 
European Space Agency (ESA) with instruments provided by two sci- culiar, spectral properties of sources in the ERCSC, which 
entific consortia funded by ESA member states (in particular the lead include peaked- spectrum, flat-spectrum, upturn- spectrum and 
countries Rrance and Italy), with contributions from NASA (USA) and multicomponent- spectrum sources. We combine archival data as 
telescope reflectors provided by a collaboration between ESA and a sci- ^^^^ ^^^^ ground-based, radio observations with the Planck 
entinc consortium led and lunded by Denmark. 
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data to construct SEDs from ~3 to -200 GHz (the exact fre- 
quency coverage varies case by case), to validate the astrome- 
try and photometry of radio sources in the ERCSC. The Planck 
data are also valuable for studying, for instance, high frequency 
peaked- spectrum sources that were previously underrepresented 
in radio source populations due to the lack of observations in 
the sub-millimeter regime. On the other hand, the ground-based 
data complement the Planck data, and are crucial in defining the 
spectral shape of some sources by extending the observed SEDs. 
We present a sample of sources with near- simultaneous Planck 
and ground-based observations (primarily employing the VLA, 
Eff'elsberg, IRAM and Metsahovi telescopes in the northern 
hemisphere, and ATCA in the southern hemisphere), to control 
for variability. We also investigate a small number of ERCSC 
sources at 30 to 70 GHz without clear identification in existing 
radio surveys. 

1.1. The Planck mission 

Planck (Tauber et al. 2010; Planck Collaboration 2011a) is the 
third generation space mission to measure the anisotropy of the 
cosmic microwave background (CMB). It observes the sky in 
nine frequency bands covering 30 to 857 GHz with high sensi- 
tivity and angular resolution from 32' to 5'. The Low Frequency 
Instrument (LEI; Mandolesi et al. 2010; Bersanelli et al. 2010; 
Mennella et al. 2011) covers the 30, 44 and 70 GHz bands with 
amplifiers cooled to 20 K. The High Frequency Instrument (HEI; 
Lamarre et al. 2010; Planck HEI Core Team 2011a) covers the 
100, 143, 217, 353, 545 and 857 GHz bands with bolometers 
cooled to 0.1 K. Polarisation is measured in all but the highest 
two bands (Leahy et al. 2010; Rosset et al. 2010). A combina- 
tion of radiative cooling and three mechanical coolers produces 
the temperatures needed for the detectors and optics (Planck 
Collaboration 2011b). Two Data Processing Centers (DPCs) 
check and calibrate the data and make maps of the sky (Planck 
HEI Core Team 201 lb; Zacchei et al. 201 1). PlancJ^s sensitivity, 
angular resolution and frequency coverage make it a powerful 
instrument for Galactic and extragalactic astrophysics as well as 
cosmology. 

The scan strategy employed in the Planck mission is de- 
scribed in Planck Collaboration (2011a). As the satellite spins, 
sources are swept over the focal plane, as indicated schemati- 
cally in Figure 1 . In the course of each day, the pointing axis of 
the telescope is adjusted by -1°, so a given source will follow a 
slightly diff'erent track across the focal plane; thus its flux den- 
sity is the average of many such scans. In addition, since the data 
included in the ERCSC amount to -1.6 full-sky surveys, some 
sources have been covered twice with a time separation of ~6 
months. Finally, sources near the ecliptic poles, where the scan 
circles intersect, are often covered multiple times. It is thus im- 
portant to keep in mind that the flux densities cited in this paper 
(and indeed in the ERCSC as a whole) are averaged. Figure 1 il- 
lustrates that flux measurements at 44 GHz are particularly sus- 
ceptible to time-dependent efl'ects, because of the wide spacing 
of the 44 GHz horns in the focal plane (see further discussion in 
§5.2.2). 



1.2. The ERCSC 

The Planck Early Release Compact Source Catalog (Planck 
Collaboration 2011c) provides lists of positions and flux den- 
sities of compact sources at each of the nine Planck frequen- 
cies. For frequencies from 30 to 143 GHz (those mostly cited 
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Fig. 1. Focal plane, showing spacing of the Planck receivers. 
Each day, the pointing is adjusted by -1° in the vertical direc- 
tion. Note that the wide separation of the three 44 GHz horns 
(two at the top, one at the bottom) causes the 44 GHz observa- 
tions of a given source to take place at two times separated by 
7-10 days for each scan. 



in this paper), sources were detected using Powell Snakes tech- 
niques (Carvalho et al. 2009). In the four highest frequency chan- 
nels, sources were located using the SExtractor method (Bertin 
& Arnouts 1996). A set of selection criteria was further applied 
to select sources that are included in the ERCSC. The primary 
criterion was a Monte Carlo assessment designed to ensure that 
> 90 % of the sources in the catalogue are reliable and have a flux 
density accuracy better than 30 %. External validation (also dis- 
cussed in Planck Collaboration 2011c) shows that ERCSC met 
its reliability criterion, and we show evidence in § 4 that the flux 
density scale of ERCSC is accurate. Secondary selection crite- 
ria, including the elimination of extended sources, were also ap- 
plied. Virtually all of the sources discussed in this paper, and the 
vast majority of extragalactic ERCSC sources, were unresolved 
by Planck and in many cases even at the much higher angular 
resolution of the VLA or other ground-based instruments. 

The flux densities in the ERCSC are calculated using aper- 
ture photometry. The efl'ective band centers or corresponding 
colour corrections depend to some degree on the spectrum of 
the source being observed. This relatively small dependence is 
discussed in the LEI and HEI instrument papers (Mennella et al. 
201 1 ; Planck HEI Core Team 201 la; Zacchei et al. 201 1 ; Planck 
HEI Core Team 2011b). We adopt here (and Hst in Table 1) 
the efl'ective central frequencies defined in those papers, and 
the colour corrections defined for a source with a spectral index 
a = -0.5 (using the convention Sy oc y^). To obtain the correct 
flux density for an assumed narrow band measurement, we di- 
vide the tabulated ERCSC flux densities by the colour correction 
factor at the corresponding central frequency. After the release 
of the ERCSC, it was found that the correction for aberration 
introduced by the motion of the satellite had not been correctly 
made. This introduces small (< 0.35') errors in the catalogued 
positions of sources; these can in turn produce very small errors 
(< 1 %) in flux densities. We have not corrected the ERCSC flux 
densities for this small efl'ect. 

In the next section, we describe briefly some general prop- 
erties of the extragalactic radio sources in the ERCSC. We 
describe in § 3 related ground-based observations. In § 4, we 
discuss the identification of sources, positional accuracy and 
flux density scales, comparing those obtained from Planck with 
ground-based measurements. This parallels the validation work 
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Table 1. Parameters of Planck bands employed in this paper. 



Planck Central Frequency Colour Beam FWHM 



Band [GHz] Correction [arcmin] 



30 28.5 1.037 32.6 

44 44.1 1.018 27.0 

70 70.3 1.031 13.0 

100 100 0.999 9.94 

143 143 1.006 7.04 

217 217 0.993 4.66 

353 353 0.990 4.41 



discussed in Planck Collaboration (201 Ic). Variability of source 
luminosity and the issue of different resolutions are also dis- 
cussed. We present and discuss in § 5 several examples of in- 
teresting classes of extragalactic radio sources. We conclude in 
the final section and point towards further research on these radio 
sources and many others contained in the ERCSC. 

2. Radio sources in the ERCSC 

The ERCSC contains hundreds of extragalactic radio sources 
at frequencies up to 143 GHz. At frequencies > 100 GHz, the 
Planck surveys are unique. At 30 and 70 GHz, the higher sensi- 
tivity and resolution off'ered by Planck allow us to detect more 
sources in a single sky survey than in seven years of survey by 
the WMAP satellite (Gold et al. 201 1). 

A major finding from the analysis of ERCSC is that many 
bright radio sources have relatively flat {a > -0.5) radio spec- 
tra extending up to and sometimes beyond 143 GHz. The vast 
majority of the extragalactic sources detected at 100, 143 and 
217 GHz are synchrotron-dominated radio sources and not 
dusty galaxies. The statistical properties of radio sources in the 
ERCSC are discussed in much greater detail in the companion 
paper, Planck Collaboration (201 Id). 

A second finding is the absence of compelling evidence for 
any new class of extragalactic radio sources. At 30 to 70 GHz, 
more than 90 % of the extragalactic sources were reliably associ- 
ated with radio sources in other large area surveys conducted at 
8-20 GHz (Planck Collaboration 2011c). Still others have plau- 
sible identifications in lower frequency radio catalogues. We dis- 
cuss those ERCSC sources without clear identification in § 5.4 
since they may potentially contain new types of sources. 

Many of the identified radio sources are blazars, some 
of them clearly variable (see a detailed discussion in Planck 
Collaboration 201 le). In § 5, we discuss a small number of ex- 
treme or unusual radio sources as additional examples of the sci- 
entific richness the ERCSC provides for the study of extragalac- 
tic radio sources. 

The Planck data used here are drawn entirely from the 
ERCSC. These data are supplemented by ground-based observa- 
tions at frequencies below and overlapping the Planck frequency 
bands, which were generally made quasi- simultaneously with 
the Planck observations of a given source (within 7-10 days, 
typically). This step was taken to monitor and control for vari- 
ability. 

3. Ground-based observations 

In planning the Planck mission, it was recognised that the sci- 
ence yield of Planck'^ millimeter and FIR sky surveys would be 
increased if accompanying ground-based observations could be 



made. These include, but are not limited to, radio observations 
at frequencies below and overlapping Planck frequency bands, 
optical observations for source identification and both ground- 
and satellite-based X- and y-ray observations. This paper dis- 
cusses only the supporting radio observations; for a discussion 
of the approximately simultaneous X- and y-ray observations, 
see Planck Collaboration (201 le). Table 2 provides information 
on the radio observatories that are involved and produced data 
used in this paper. Most of the radio observations were made 
preemptively — that is, we observed sources expected to be de- 
tected by Planck at about the same time they passed through the 
Planck beam. The POFF software (Massardi & Burigana 2010) 
was used to predict which sources would be seen by Planck in 
a given week. Obviously, only previously known sources can be 
observed preemptively. We later conducted a small number of 
follow-up observations on ERCSC sources. 

In the southern hemisphere, a substantial amount of time was 
obtained at the Australia Telescope Compact Array (ATCA) to 
make preemptive observations at a wide range of frequencies 
up to and overlapping with the Planck. The Planck- ATCA Co- 
eval Observations (PACO) project (Massardi et al. 2011) con- 
sists of several-epoch observations of a compilation of sources 
selected from the Australia Telescope 20 GHz survey (AT20G, 
Murphy et al. 2010) with \ b\ >5° . Observations were made with 
ATCA in the frequency range between 4.5 and 40 GHz at epochs 
close in time to the Planck observations over the interval from 
July 2009 to August 2010. The PACO sample includes a com- 
plete flux density limited and spectrally selected sample over the 
whole Southern sky; 147 PACO point-like sources have at least 
one observation within 10 days of the Planck observations. 

In the northern hemisphere, quasi- simultaneous cm/mm ra- 
dio spectra for a large number of Planck blazars have been ob- 
tained within the framework of a Fermi-GST related monitor- 
ing program of y-ray blazars (F-GAMMA program, Fuhrmann 
et al. 2007; Angelakis et al. 2008). The frequency range spans 
from 2.64 GHz to 142 GHz using the Efl'elsberg 100 m and 
IRAM 30 m telescopes. The Efl'elsberg measurements were con- 
ducted with the secondary focus heterodyne receivers at 2.64, 
4.85, 8.35,10.45, 14.60, 23.05, 32.00 and 43.00 GHz. The obser- 
vations were performed quasi- simultaneously with cross-scans, 
that is slewing over the source position in azimuth and elevation 
direction with the number of sub- scans adjusted to reach the de- 
sired sensitivity (for details, see Fuhrmann et al. 2008; Angelakis 
et al. 2008). Subsequently, pointing ofl'-set correction, gain cor- 
rection, atmospheric opacity correction and sensitivity correc- 
tion were applied to the data. The IRAM 30 m observations were 
carried out with calibrated cross-scans using the EMIR horizon- 
tal and vertical polarisation receivers operating at 86.2 and 142.3 
GHz. The opacity corrected intensities were converted into the 
standard temperature scale and finally corrected for small re- 
maining pointing off'sets and systematic gain-elevation eff'ects. 
The conversion to the standard flux density scale was done us- 
ing the instantaneous conversion factors derived from frequently 
observed primary (Mars, Uranus) and secondary (W3 (OH), K3- 
50A, NGC7027) calibrators. 

In the northern hemisphere, some data also come from the 
Metsahovi telescope operating at 37 GHz. The observations 
were made with the 13.7m Metsahovi radio telescope, which 
is a radome-enclosed, paraboloid antenna situated in Finland. 
The measurements were made with a 1 GHz-band dual-beam re- 
ceiver centered at 36.8 GHz. The observations are ON-ON, al- 
ternating the source and the sky in each feed horn. A typical inte- 
gration time to obtain one flux density data point is 1200-1400 s. 
The detection limit of the telescope at 37 GHz is of the order of 
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Table 2. Ground-based radio observations. 



Observatory 


Project Leaders 


Frequencies [GHz] 


ATCA, Australia 


Massardi 


4.5-40 


Effelsberg, 


Fuhrmann, Angelakis 


2.6-43 


Germany 


and Rachen 




IRAM, Spain 


Fuhrmann, 


86.2, 142.3 




Ungerechts and 






Rachen 




Metsahovi, 


Lahteenmaki and 


37 


Finland 


Tornikoski 




VLA/EVLA, 


Partridge and Sajina 


5, 8, 22, 43 


USA 







0.2 Jy under optimal conditions. Data points with a signal-to- 
noise ratio < 4 are treated as non-detections. The flux density 
scale is set by observations of DR21. Sources NGC7027, 3C274 
and 3C84 are used as secondary calibrators. A detailed descrip- 
tion of the data reduction and analysis is given in (Teraesranta 
et al. 1998). The error estimate in the flux density includes the 
contribution from the measurement rms and the uncertainty of 
the absolute calibration. 

In addition, small amounts of time, scattered throughout 
the first 17 months of the Planck mission, were obtained at 
the Very Large Array (VLA) of the National Radio Astronomy 
Observatory (NRAO). Observations at the VLA began 24 July 
2009, slightly before the beginning of Planck'^ first sky sur- 
vey. The first set of observations on 24 July 2009 is the only 
set not approximately simultaneous with Planck observations. 
VLA/EVLA measurements continued at irregular intervals un- 
til November 2010, with a substantial gap in the spring of 2010 
when the VLA was converted to the EVLA. Most of the VLA 
and EVLA runs were brief 1-2 hour chunks of time; 5-8 Planck 
sources were typically observed per hour, besides flux calibra- 
tors and phase calibrators. In many cases, VLA flux and phase 
calibrators were of interest themselves because they were bright 
enough to be detected by Planck. The integration times were 
typically 30 seconds at 4.86 GHz and 8.46 GHz, 100 seconds at 
22.46 GHz, and 120 seconds at 43.34 GHz. AH u-v data were 
flagged, calibrated and imaged using standard NRAO software: 
AIPS or CASA. The flux density measurements were calibrated 
against one or both of the primary calibrators used by NRAO: 
3C48 or 3C286. 

The positions and flux densities of sources observed by the 
VLA are listed in Table 3. Since the VLA flux densities are accu- 
rate to a few percent for these bright sources, individual uncer- 
tainties are not included. Typical 1 a errors at 22 GHz ranged 
from 2 mJy (for the fainter sources) to 15 mJy (for brighter 
sources); and from 3-15 mJy at 43 GHz; both are small com- 
pared to the uncertainties in Planck flux densities. The VLA 
(EVLA) was in difl'erent configurations at diff'erent times; hence 
the angular resolution of the array was changing. Additionally, 
for a given configuration, the resolution was much finer at higher 
frequencies. We thus flagged sources that showed signs of res- 
olution in any configuration at any frequency. In general, the 
VLA observations were timed to occur during the same week 
as Planck was expected to see a given source, but that was not 
always possible. Therefore, each observation is tagged with a 
date in Table 3, to be compared with the dates of observation 
tabulated in the ERCSC. Given that the typical timescale for 
variability of radio sources is roughly weeks to months (see 



e.g., Hovatta et al. 2008; Nieppola et al. 2009), we were able to 
control for variability to some degree. Nevertheless, it must be 
borne in mind that the ground-based observations were not al- 
ways taken at the same time Planck was observing a source. The 
Planck beams at diff'erent frequencies also swept over a source at 
diff'erent times due to the extended layout of its focal plane. We 
use this set of data in § 4 to check the positional and photometric 
accuracy of ERCSC sources. 



4. Identification and validation 

To construct SEDs that include measurements from both Planck 
and ground-based telescopes, we need to ensure that the identi- 
fication of a source in the ERCSC has been correctly made, and 
that the flux density scales are consistent. Both properties were 
extensively checked as part of the validation work of the ERCSC 
(Planck Collaboration 2011c). We present here the comparison 
of quasi- simultaneous VLA and Planck observations of a set of 
bright radio sources. 

4.1. Source Identification 

Identifying a source for supporting ground-based observations 
requires an understanding of the astrometric accuracy of the 
ERCSC. In Figure 2, we show a comparison of ERCSC source 
positions at 30 and 70 GHz with those measured for the pre- 
sumed identification by the VLA at 22 GHz. The VLA positions 
for such bright sources are typically accurate to a few arcsec- 
onds even for the most compact configurations. Figure 2 indi- 
cates that (a) the ERCSC meets its specification of having a po- 
sitional accuracy good to FWHM/5 (or -6.5' at 30 GHz); (b) a 
search radius of 1/2 FWHM of the beam around each ERCSC 
source position is sufficient to locate any related source; and 
(c) any bright radio source within a few arcmin of the ERCSC 
position is likely the correct low frequency counterpart. This 
last point can be made more quantitative using the AT20G sur- 
vey, which covers the entire southern hemisphere to ~40 mJy at 
20 GHz, and is essentially complete above 100 mJy. The density 
of sources with S > 100 mJy is -0.25 per square degree, and the 
integral counts of sources have a slope of -1.15. Thus we ex- 
pect on average 0.07 sources-deg"^ above 300 mJy, the thresh- 
old we applied when cross-correlating AT20G sources with the 
ERCSC. At 30 GHz, the probability of a random AT20G source 
appearing within 16.25' of a given ERCSC source is ~ 1 .6 %. The 
ERCSC contains several hundred extragalactic sources, depend- 
ing on frequency; only 5-8 of these might be falsely identified 
or confused. Additionally, all the sources described in § 5 fall 
within a radius of 5' of the ERCSC position, further reducing the 
chance coincidence rate by a factor of -10. The ERCSC source 
positions were also independently compared to those recorded 
for several hundred bright quasars. The median scatters in off'set 
are 2.0', 1.7', 1.1', 0.8', 0.7' and 0.3' for the Planck frequency 
bands 30 to 217 GHz (Planck Collaboration 201 Ic). This is con- 
sistent with our observations shown in Figure 2. 

4.2. Comparison of flux density scales 

The flux density scale of Planck is ultimately tied to the ampli- 
tude of the CMB dipole, and careful measurements of the an- 
gular resolution of Planck and the ground-based instruments are 
required in order to convert from temperature units to flux den- 
sity (Planck Collaboration 201 Ic). Since the calibration standard 
was the very large angular scale dipole signal, it is important to 
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Table 3. VLA observations on a set of bright ERCSC sources. 



ERCSC name 


VLA RA 


VLA Dec 


S 4.86 GHz 


S 8.46 GHz 


S 22.46 GHz 


S 43.34 GHz 


VLA obs. date 


Planck obs. date'' 




[J2000] 


[J2000] 


[Jy] 


[Jy] 


[Jy] 


[Jy] 


[yymmdd] 


[yymmdd] 


G052.38-36.49 


21:34:10.3 


-01:53:17 






1.698 


1.218 


091103 


091102, 091103, 091108-091110, 


















100509-100511, 100517 


0056.70+80.65^ 


13:31:08.3 


30:30:33 


7.485 


5.202 


2.516 


1.444 


091022, 100103 


091230, 091231 


G075.68-29.62 


22:03:26.9 


17:25:48 


0.987 


1.055 


1.035 


1.073 


091211 


091121, 091128-091202, 100524-100527, 


















100602, 100603 


G077.45-38.54 


22:32:36.4 


11:43:51 


5.594 


5.827 


4.407 


3.606 


091211 


091127, 091203-091207, 100531-100603 


















100506-100510, 100519, 100520 


G085.72+26.08 


18:24:07.1 


56:51:01 


1.403 


1.876 


1.771 


1.258 


090724 


090827-090829, 100417-100422, 


















100503,100504 


G085.86+83.31 


13:10:28.7 


32:20:44 


1.427 


2.133 


2.856 


1.977 


100103 


091221-091225, 100102, 100103 


G085.95-18.77^ 


22:03:15.0 


31:45:38 


2.165 


2.495 


2.576 


2.841 


091211 


091202, 091203, 091211-091215, 


















100531-100603 


G086. 12-38. 18^ 


22:53:57.8 


16:08:54 


9.463 


8.176 


11.834 


21.570 


091211 


091206, 091207, 091213-091216 


G090. 11-25.64 


22:36:22.5 


28:28:57 


1.285 


1.273 


1.412 


1.547 


091211 


091210, 091211, 091218-091221 


G090. 14+09.66 


20:23:55.8 


54:27:36 


0.952 


1.132 


0.804 


0.502 


090724 


091222, 091223, 100517, 100518 


G092.62-10.44 


22:02:43.3 


42:16:39 


3.994 


3.996 


3.440 


3.085 


091211 


091213, 091214, 091222-091226, 100606, 


















100607 


G096.08+ 13.77 


20:22:06.7 


61:36:59 


3.110 


3.800 


1.946 


0.836 


090724 


091220-091222, 100103-100108, 


















100526-100531, 100602 


G097.46+25.04 


18:49:16.1 


67:05:42 


1.284 


2.431 


3.553 


2.956 


090724 


091012-091014, 091017-091022, 


















091112-091114, 091225-091227, 


















100112-100119, 100512-100514, 


















100516-100521 


G098.28+58.29 


14:19:46.6 


54:23:14 


1.178 


1.249 


1.231 


1.176 


100103 


091206, 091207, 091209-091212, 091226, 


















091227 


G098.51+25.79 


18:42:33.6 


68:09:25 


0.765 


0.892 


0.676 


0.461 


090724 


091013-091015, 091110-091112, 091231, 


















100101, 100116-100119, 100121, 100122, 


















100514-100516, 100518-100521, 100523, 


















100524 


Gl 10.05+29.07 


18:00:45.7 


78:28:04 


2.098 


2.850 


2.598 


1.936 


090724 


091011-091016, 091026, 091027, 100202, 


















100203, 100214-100219 


G126.44-64.29" 


00:57:34.9 


-01:23:28 


0.338 


0.097 


0.060 


0.070 


091103 


100107, 100108 


G128.95+11.97 


02:17:30.8 


73:49:33 


4.168 


4.090 


2.623 


1.709 


090827 


090911-090914 ,090920, 090921, 100211, 


















100212, 100217-100220 


G133.94-28.63^ 


01:37:41.3 


33:09:35 


5.426 


3.140 


1.098 


0.501 


090827, 091211 


100130, 100131 


G143.55+34.41 


08:41:24.4 


70:53:42 


1.762 


1.485 


2.267 


2.900 


091022 


091012-091014, 091019, 100316, 100317, 


















100324-100327 


G145.58+64.96 


11:53:24.5 


49:31:09 


1.167 


0.949 


0.885 


0.798 


100103 


091116-091119, 091125, 091126, 100511, 


















100512, 100522-100526 


G174.47+69.79 


11:30:53.3 


38:15:19 


1.334 


1.356 


1.019 


0.680 


100103 


091121, 100527 


G195.26-33.13 


04:23:15.8 


-01:20:33 






7.500 


7.289 


091103 


090828-090830, 090903, 100221, 


















100225-100227 


G199.42+78.39 


11:59:31.8 


29:14:44 


1.609 


1.530 


1.416 


1.141 


100103 


091204-091207, 091214, 100531, 100601 


G200.06+31.89 


08:30:52.1 


24:11:00 


1.269 


1.246 


1.050 


0.825 


091022 


091021-091023, 091027, 100409, 100410, 


















100414-100416 


G208.18+18.75 


07:50:52.0 


12:31:05 


3.685 


4.486 


4.400 


3.939 


091022 


091014-091016, 091020, 100404, 


















100408-100410 


G211.33+19.05 


07:57:06.6 


09:56:35 


1.143 


1.530 


1.863 


1.846 


091022 


091015-091018, 091022, 100406, 


















100410-100412 


G216.97+11.36 


07:39:18.0 


01:37:04 


1.143 


1.260 


1.514 


1.553 


091022 


091012-091014, 100408-100410 


G221.26+22.36 


08:25:50.3 


03:09:24 


0.777 


0.737 


0.668 


0.648 


091022 


091024, 091025, 100419, 100420 


G251.59+52.70 


10:58:29.6 


01:33:59 


3.193 


4.186 


5.317 


5.078 


100103 


091207-091210, 091216, 100524, 


















100530-100602 


G255.00+81.65 


12:24:54.5 


21:22:47 


1.233 


1.267 


1.412 


1.184 


100103 


091218-091221 


G283.75+74.54^ 


12:30:49.4 


12:23:28 


62.000 


39.000 


12.000 




100103 


091226-091229, 100104, 100628, 100629 



Notes. Planck observation dates record the days that the source is observed by any of the 30, 44 and 70 GHz channels. Standard VLA flux 
calibrator, observed twice. Reported flux densities are the average of the two observations. Resolved source in the VLA. 



confirm that flux density measurements of compact, unresolved, 
sources are accurate. 

A potential complication in comparing flux densities is the 
greatly diflTerent angular resolution of Planck and the ground- 
based instruments. As discussed above, depending on the fre- 
quency and configuration, the angular resolution of the VLA can 
be as small as a fraction of an arcsecond, whereas the Planck 
beams are typically many arcminutes in size. Thus the VLA 
can "miss" flux included in the Planck beam. To first order, we 
take account of this by excluding from the flux density com- 
parison any source seen to be resolved at any frequency by the 



VLA (see Table 3). A clear example of such a source is 3C274 
(J1230-hl223). 

We then compare in Figure 3 the VLA measurements with 
the ERCSC flux densities at 30, 44 and 70 GHz. The VLA 22 
and 43 GHz measurements were used to interpolate the VLA 
flux densities to the center frequency of Planck 30 GHz band. 
In the case of 44 GHz, the VLA measurements centered at 43.34 
GHz and the Planck measurements centered at 44.1 GHz were 
directly compared without any correction. We performed least 
squares fits to the data in the linear regime. The agreement at 
30 GHz (28.5 GHz to be precise) is very good, with the mea- 
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70GHz 
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Offset in RA [orcmin] Offset In RA [orcmin] 

Fig. 2. Comparison of the VLA measurements made at ~22 GHz with the ERCSC source positions at 30 GHz (left) and 70 GHz 
(right). 





Fig. 3. Upper left: Comparison of Planck measurements in the "30 GHz" band with interpolated values from the VLA. Upper right: 
Planck vs. VLA measurements at 44 GHz. Lower center: Planck measurements at 70 GHz compared to the extrapolated VLA 
measurements. Note that although the log-log plot is presented for clarity, the fits were done on the linear data. 



sured slope equal to LOS and an intercept of 0.08. At 44 GHz, 
the agreement is also reasonably good, with a measured slope of 
LOl and an intercept of 0.24. 

At 70 GHz, we extrapolate the VLA measurements from 43 
GHz using the 22-43 GHz spectral index. As shown in Figure 3, 
the extrapolated VLA values diff'er from the Planck values by 
~20%. There is strong evidence, presented in the companion 



paper (Planck Collaboration 201 Id), that the spectral index of 
radio sources detected by Planck steepens at frequencies above 
70 GHz or perhaps 44 GHz for some sources. Thus some of our 
extrapolated VLA 70 GHz measurements may be biased high. 
Green Bank Telescope (GBT) 90 GHz measurements of a sam- 
ple of VLA sources (Sajina et al. 2011), show that spectral cur- 
vature indeed plays a role. Among our VLA sample, 2 sources 
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have 90 GHz measurements. Replacing the extrapolated 70 GHz 
values with interpolated ones from VLA 44 GHz and GBT 90 
GHz changes the slope of the best fit line to be closer to unity, 
from 0.78 to 0.86. We emphasise that the 90GHz observations 
were not made simultaneously with the Planck measurements, 
so they are less well controlled for variability than the quasi- 
simultaneous VLA (or EVLA) observations. 

We also examined the median ratios of the Planck measure- 
ments and the extrapolated VLA flux densities. These are respec- 
tively L08 ± 0.04, L12 ± 0.07 and 0.99 ± 0.06 for the 30, 44 and 
70 GHz bands. These median values have the advantage over the 
linear fits in that they eff'ectively force a zero intercept and mit- 
igate the eff'ect of outliers. The median ratios confirm the good 
agreement at 30 and 44 GHz and show that the agreement be- 
tween our VLA and Planck flux densities at 70 GHz is better 
than implied by the linear fit. There is still, however, substantial 
scatter with the standard deviations of these ratios being 0.22, 
0.38 and 0.29 for the 30, 44 and 70 GHz bands, respectively. 
The scatter is likely dominated by the intrinsic variability of our 
sources. 



J2022+6137 



10.0 
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V [GHz] 

Fig. 4. SED of a known CSO source, J2022-h6137. Planck mea- 
surements are shown in red, WMAP-lyr flux densities in blue, 
our new VLA measurements in pink and archival data obtained 
from NED in green. 



5. Extreme radio sources 

In this section, we present several examples of sources illustrat- 
ing the broad range of high frequency spectral behavior in bright, 
extragalactic radio sources as observed by Planck This is not in- 
tended as an exhaustive list of extreme sources or an extensive 
discussion to explore the phenomenology of these sources, but 
rather a sample from the rich data set provided in the ERCSC. 
Here we employ both Planck and ground-based observations, 
using the latter to extend the source SEDs to lower frequen- 
cies. In some cases, the ground-based observations were made 
at roughly the same time as Planck observed the source, as men- 
tioned in § 3. Archival data were used when such nearly simul- 
taneous observations were not available. 

5.1. Peaked spectrum sources 

The term Gigahertz-Peaked- Spectrum (GPS) sources, in princi- 
ple, refers to a morphological type of the radio spectra, namely, 
to a spectral index > for y < Vp, and a < for y > Vp with 
a peak frequency Vp in the GHz regime. The GPS phenomenon 
was originally thought to be associated with compact, putatively 
very young radio sources (O'Dea 1998), and in fact many exam- 
ples for this association have been found (Conway 2002). Recent 
research, however, has shown that a large fraction of sources 
with GPS features are not of this type, but associated with com- 
pact, beamed jet sources, commonly identified as blazars. These 
two classes produce slightly diff'erent spectra: the former shows 
a narrow peak, while the latter is typically broadly peaked. To 
date, however, the only secure method to distinguish between 
them is to employ VLBI observations, revealing the source mor- 
phology (Bolton et al. 2006; Vollmer et al. 2008). In this paper, 
we use the term "GPS sources" purely phenomenologically, and 
discuss in the following sections the diff'erent, known, classes of 
sources for which we have examples in the ERCSC and also a 
set of ERCSC sources that are candidate GPS sources. 

5.1 .1 . ERCSC spectra of known Compact Symmetric 
Objects 

One class of GPS sources is Compact Symmetric Objects (CSO), 
thought to be either very young or very recently activated radio 



galaxies (Owsianik & Conway 1998). Mapped with VLBI res- 
olution, these sources show a typical, symmetric radio-double 
morphology, but with linear extensions of 10 pc or less. They 
produce GPS type spectra by essentially single-zone synchrotron 
emission with synchrotron self- absorption causing a spectral 
turn-over at v > 1 GHz (O'Dea 1998). Unlike blazars, their 
emission is most likely not Doppler boosted, and they show no 
hint of fast variability, although their interpretation as young ra- 
dio galaxies suggests that a spectral evolution over time scales of 
decades may be possible. There is also a well-established con- 
nection between the peak frequency of a CSO-type GPS source 
and its linear extent: d ~ lOOVp^pc, with Vp in GHz. Thus 
sources with peak frequencies in the Planck regime would point 
to objects at most a few parsecs in linear size. Previously, pop- 
ulations of sources peaking above an observed frequency of 10 
GHz (in the observer's frame) might have been strongly under- 
represented because they are relatively faint at the low frequen- 
cies where most large surveys have been made. Planck opens 
the possibility of detecting such "extreme GPS" sources, or 
"high frequency peakers" in the nomenclature of Dallacasa et al. 
(2000). In addition, Planck allows the examination of the spec- 
tral decline of GPS sources at frequencies far above the peak fre- 
quency where optically thin synchrotron radiation is thought to 
dominate the emission. In Figure 4, we show the spectrum of one 
known CSO source, J2022-h6137 (B2021-h614 from Conway 
2002). Planck data have been combined with archival data as 
CSO sources are not expected to be variable. 

5.1 .2. Additional high frequency peakers in ERCSC 

As mentioned above, Planck observations make possible the de- 
tection and confirmation of rare, bright, GPS sources with very 
high peak frequencies. We provide in Table 4 a list of poten- 
tial GPS sources that show convex SEDs constructed using the 
30 to 143 GHz ERCSC data, along with the archival data at 20, 
8.6 and 4.8 GHz (Murphy et al. 2010; Healey et al. 2007) if 
available. We restricted our search to | | > 10° sources. While 
all of the tabulated sources display spectral peaks in the Planck 
frequency range, there is no guarantee that they are particularly 
young and compact CSO sources. Indeed, some of these bright 
sources have been previously studied, and most of them appear 
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Fig. 5. SEDs of sources that show spectral peak in the Planck bands. ERCSC data are shown in filled circles and low frequency 
archival data are shown in open circles. 
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Fig. 6. With the additional Planck data, source J1604+5714 is a 
newly disclosed GPS source candidate that previously showed 
an inverted spectrum at low frequencies. ERCSC data are shown 
in filled circles and low frequency archival data are shown in 
open circles. 

to be variable flat spectrum sources (see the notes in Table 4 and 
Torniainen et al. 2005). We discuss sources of this kind in more 
detail in § 5.2. 

Figure 5 shows four sources that are quite undistinguished 
in low frequency catalogues, but reveal peaks in the Planck fre- 
quency range. Note that sources J0639+7324, J0754+4820 and 
J2 100-2932 all present peaks at 70 GHz or higher. The recorded 
redshift of J0639+7324 is 1.85; thus the rest frame frequency of 
the emission peak is extremely high at ~200 GHz. In addition. 



Planck data also helped to disclose the GPS-like spectrum for 
some sources that show inverted spectra at lower frequencies, as 
in the case of J1604+5714 (Figure 6). 

As already mentioned, it is not possible from spectral infor- 
mation alone to identify high frequency peakers in the ERCSC 
with a new, very compact, radio source population. In fact, all 
but 5 sources in Table 4 can be identified with known blazars 
in the Roma-BZCAT blazar catalogue (Massaro et al. 2008). 
Out of the five exceptions, 3 are further associated with blazars 
in the CGRaBS blazar catalogue (Healey et al. 2008). The last 
two sources are identified as NGC1218 (J0308-h0405) and HB89 
2002-185 (J2005-1821)in NED. The SEDs of these two sources 
are shown in Figure 7. NGC1218 has a typical steep spectrum of 
radio galaxy up to 20 GHz, then a bump is seen at -40 GHz. 
Since WMAP (in the 7 year co-added map) did not catch this 
feature whereas Planck did (during its 1 .6 sky surveys), this peak 
is clearly caused by variability, suggesting some blazar like fea- 
tures in the galaxy, likely from close to the center. The spectral 
feature of HE 89 2002-185 is evidently similar to that of source 
J2 100-2932 in Figure 5 which has been identified as a flaring 
blazar. Therefore we suggest that this source could very well be 
a blazar. Further observations are needed to confirm this assump- 
tion. Our finding that the peaked spectra in many bright sources 
is the result of flaring is consistent with earlier discussions by 
Tornikoski et al. (2001); Torniainen et al. (2005); Bolton et al. 
(2006). 
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Table 4. GPS source candidates in the ERCSC. 
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1 OT 1 oc 

i.U/+U.Uj 


f\ A A A 000 

U.444UUU 


B 


JU4jj— 4DiD 


rv/i .cccQ rv 
U4. JJ. JJ.U 


— 4d. io. JU 


1 rvOj^rv rv? 
i.UZ+U.U/ 


1 11 j^rv rvrv 

i.ji+u.uy 


1 'TOj^rv 1 c 
i. /y+U.i J 


01 j^o 1 
Z.Ui+U.iZ 


^Qj^O 0^ 

Z.dJ+U.Ud 


/I 1 j^o 1 
4.iD+U.Zi 


Q ^1 j^o 1 
J.oi+U.io 


^1 j^O 1 Q 

Z.oi+U.i J 


ocoooo 
U.ojZoUU 


B 


JU4j /+Ud4U 


rv/i . c?. 1 /I ^ 
U4. J /. i4.D 


UD.4U.4i 


rv ccj^rv rvc 
U.jj+U.Uj 


rv OQj^rv rvo 
U.oJ+U.Uo 


1 rvrvj^rv i o 
i.Uy+U.io 




OOj^O 1 

U.oU+U. iU 




/I Q 
U.4J 


^0 
U.oZ 


/lOCOOO 

U.4UjUUU 


B 


JUjZj— Zjjo 


rvcocrv/i o 
Uj.Zj.U4.o 


OQ.Q^./irv 
— ZJ. JD.4U 


rv /lOj^rv rvQ 
U.4Z+U.UJ 


rv ^Cj^rv rvT 
U.OJ+U.U/ 




1 Q'Tj^O 1 ^ 

i.J /+U.iD 


1 1 Qj^O 0^ 

i. i J+U. Uo 


'toj^o oc 
U. /y+U.Uj 


ooj^o oc 
U.oy+U.Uj 


001 j^O O/I 

U.oi+U.U4 


Q 1 00000 
J.iUUUUU 


B 


JUjZd— 4oJi 


rvco/c.rvrv n 


/I o.Q 1 .rvo 
— 4o. Ji.Uo 


rv QOj^rv rv/c 
U.jy+U.Uo 


rv ^Cj^rv rvT 
U.OJ+U.U/ 




1 O^j^O 1 

i.yo+u.iz 


1 QOj^O 0^ 

i.JU+U.Uo 


ooj^o 01 

u.zy+u.ui 


Q/i j^o 00 
U.J4+U.UZ 


qtj^o 00 
U.J /+U.UZ 




B 


JUjj4— DiU/ 


Uj. j4. JO. / 


^1 .rvv.o^ 
— Oi.U/.Zo 




rv coj^rv rvQ 
U.jZ+U.UJ 




OQ j^O 1 

U.oJ+U.iU 


TOj^O O/I 

U. /9+U.U4 


r\ A A j^o 00 
U.44+U. UZ 


/I c j^o 00 
U.4J+U.UZ 


r\ A A j^o 00 
U.44+U. UZ 


1 OOTOOO 

i.yy /UUU 


B 


JUjJo-44Uj 


Uj:3o:j4. / 


A A .rvc . 1 Q 

-44:Uj:ij 


1 rv rvT 
D.D/+U.U/ 


/I c 1 rv rv^ 
0.4j+U.Ud 


rv 10 1 rv 1 'J 
9.JZ+U.1 J 


/I 1 1 1 
0.4D+U.1 J 


00 1 O'T 

9.UU+U.U/ 


c 00 1 oc 
J.Z9+U.ZJ 


A OQ 1 1 

4.ZJ+U.Zi 


Q 00 1 1 
J.oU+U.i9 


OO/IOOO 

U.094UUU 


TD \T 

B, V 


JUj4U— j4io 


rvc/irv.C7 i 
UJ.4U.J /.i 


c/i . 1 0.00 
— j4. io.Zy 


rv /I Q_i_rv rvQ 
U.4J+U.UJ 


rv 'TO_i_rv rvc 
U. /Z+U.Uj 


rv o'T_i_rv 1 rv 
u.y /+U.iU 


1 1 n_i_n no 
i.iU+U.Uo 


1 o^_i_o oc 
i.Uo+U.Uj 


1 1 Q_i_o 0^ 
i.i J+U. Uo 


oo_i_o oc 

u.yy+u.uj 


T 1 _i_0 O/I 

U. /i+U.U4 


1 1 ocooo 
i.iojUUU 


B 


JUjjU— J /JZ 


rvccrv.rvT a 
Uj.jU.U/.U 


CT.QO.rv^ 


rv ^rv_i_rv rvo 


f\ n A _i_rv rv/i 
U. /4+U.U4 


1 rvi _i_rv rvT 
i.Ui+U.U/ 


1 C/I _i_o 00 
i.j4+U.Uo 


1 A C_i_0 OC 

i.4j+U. Uj 


1 oo_i_o oc 
i.UU+U.Uj 


1 oQ_i_o oc 
i.UJ+U.Uj 


OQ_i_o oc 
U.yJ+U.Uj 


001 000 


B 


JUoZZ— d4Jj 


rv^.oo.cc 
Uo.ZZ.jj.y 


— 04. Jj. i / 


rv CQ_i_rv rvQ 
U.jJ+U.UJ 


rv 'T^_i_rv rvQ 
U. /D+U.UJ 


rv o^_i_rv rv^ 
U.OD+U.UO 


1 0'T_i_0 O'T 

i.U/+U.U/ 


00_i_0 OQ 

u.^^y+u. UJ 


'T0_i_O O/I 

U. /Z+U.U4 


O/I _i_0 O/I 

U.o4+U.U4 


0T_i_O O/I 

U.o /+U.U4 


1 00000 
U.iZoooy 


B 


JUdJ'5'+ /JZ4 


rv/^.QO. 1 rv 1 
UoiJyiiU.i 


.o/i .orv 

/J:Z4:Zb' 


rv rvrv i rv rv/i 
U.9U+U.U4 


1 O'T 1 rv rv/^ 
1.Z/+U.UD 


1 irv 1 rv 1 'J 
1.J9+U.1 J 


1 1 /I 1 1 

1.14+U.19 


001 1 O/^ 

U.oi+U.Uo 






U. / J 


T 1 
U. /i 


1 ocoooo 
I.ojUUUU 


B 


JUD4D+44ji 


rv^./i/c.Qi ^7 
Ud.4d. Ji. / 


/I /I . c 1 .CO 

44. ji. JO 


rv rvOj^rv rv/c 
U.yo+U.Uo 


1 o'Tj^rv rvrv 

i.z/+u.uy 


1 o^j^rv 1 rv 
i.oD+U.i9 


1 'TCj^O 1 

i. /j+U.io 


/coj^o 00 
Z.oZ+U.Uo 




0^ 
Z.Zo 


1 00 

i.ZZ 


Q QO/COOO 

J.JyoUUU 


T3 \7 
B, V 


JU/Ui— 4dJJ 


rvT ^^\^ • A A o 
U/.Ui.44.Z 


— 40. JJ.Zy 


rv ^^_i_rv rvQ 
U.DD+U.UJ 


rv o^_i_rv rv/i 
U.OD+U.U4 


1 rvo_i_rv rvo 

i.uz+u.uy 


1 Q/i _i_n no 
i.J4+U.Uo 


1 oc_i_o oc 
i.Uj+U.Uj 


1 0'T_i_0 OC 

i.U/+U.Uj 


0^_i_0 O/I 

U.OD+U.U4 


CC_i_0 OQ 

U.JJ+U.UJ 


000000 
U.oZZUUU 


B 


JU/i /+4jjy 


rv?. 1 T. crv 
U/. i /. jU.Z 


A CQO./I'T 

4j.jy.4/ 


rv QOj^rv rv/c 
U.Jo+U.Uo 




n ooj^o 1 /I 
U.^o+U. 14 




oCj^o 00 
U.yj+U.Uo 




C/^ 

U.JO 


/IT 

U.4/ 


O/I 0000 
U.y4UUUU 


B 


JU/jU+iZJi 


rvT.crv.crv a 
U/. jU. jU.4 


1 O.Q 1 .Q/l 

iZ. Ji. J4 


rv/i j^rv rvQ 
Z.U4+U.UJ 


^rvj^rv rvT 
Z.oy+U.U/ 


Q 1 Oj^O 1 c 

J.io+U.i J 


Q OOj^O 1 c 

j.yu+u.i J 


/I C/I j^o 00 
4.j4+U.Uo 




1 OT 


1 O/I 

i.Z4 


000000 
U.ooyUUU 


T3 \T 

B, V 


JU/j4+4oZU 


rvT . c /I . /1 rv 

U/:j4:4o.U 


A o.orv. 1 T 

4o:ZU:i / 


rv /^c 1 rv rvT 
U.DJ+U.U/ 


rv rv/i 1 rv rvo 
U.94+U.U0 


rv orv 1 1 'T 
U.9U+U.1 / 




1 00 
U.09+U.U9 




00 

u.zu 


01 
U.ZJ 


ITT 1 /I 

U.J / /i4Z 


B 


J Uo JO— ZZ J / 


rvo.Q/c.co Q 
Uo. JO. jZ.J 


00. Q'?. rvi 
— ZZ. J /.Ui 


rv QOj^rv rv/c 
U.jy+U.Uo 


rv /1 j^rv rvo 
U.4o+U.Uo 


n ooj^o 1 
U.oU+U. iZ 




1 /IOj^O O'T 

i.4U+U.U/ 


/l^j^O OQ 

U.4D+U.UJ 


qcj^o 00 
U.Jj+U.UZ 


Q 1 j^o 01 
U.Ji+U.Ui 


OQTOOO 

U.oJ /UUU 


B 


inO/l 1 1 'TACQ 

JUo4i+ /UjJ 


rvo. /1 1.100 
Uo.4i. iZ.Z 


'vrv. CQ .O/I 
/U. jJ.Z4 


1 rvj^rv rv/i 
Z.iU+U.U4 


rvi j^rv rvT 

z.yi+u.u/ 


Q c^j^n 1 /I 
J.JD+U.i4 


Q 1 'Tj^O 1 

J.i /+U.iZ 


OOj^O 0^ 

Z.oZ+U.Uo 




1 TC 
i. /j 


Q/i 
Z.J4 


1 TOOOO 

Z.i /ZUUU 


T) 

B 


JUo4/— Uojy 


rvo. /I'?. CO ^ 
U0.4/.JZ.D 


rv^. crv.rvo 
— Uo. jy.UZ 


rv c/i j^rv rvc 
U.j4+U.Uj 


rv rvOj^rv rv^ 
U.yZ+U.Uo 


n ooj^o 1 /I 
U.yy+U. 14 




001 j^O 0^ 
U.yi+U.Uo 


'TQj^O O/I 

U. /J+U.U4 


/IT 

U.4/ 


/I /I 
U.44 




B 


JUoj4+ZUUd 


rvo . c /I . /I Q rv 

Uo:j4:43.b' 


orv.rv/^ .rv/i 

ZU:Uo:U4 


c /irv 1 rv rvrv 
J.4U+U.U9 


^ o/^ 1 rv rvo 
D.Zo+U.Uo 


/IO 1 1 

D.49+U.1D 


'T Q /I 1 1 

/.J4+U.1D 


'T 1 1 00 

/.io+U.Uo 




/1 1 
J.4i 


ni 


OO/^OOO 

U.JUoUUU 


B 


JUyZU+444i 


rvrv. orv. CO o 
Uy.ZU. JO.o 


A A . A^ .^ 1 
44.4i. i J 


1 /I Q j^rv rv/c 
i.4J+U.UD 


1 QOj^rv rvT 

i.jy+u.u/ 


1 01 j^O 1 c 
i.oi+U.i J 


1 Cj^O 1 c 

z.i J+U.i J 


ooj^o 00 
Z.Zy+U.Uo 




1 'JT 
i.J / 


1 00 

i.U9 


1 0001 
Z.ioyyiU 


B 


JU'5'ZJ+Zol / 


rvrv . . c /I rv 

Ub':ZJ:j4.U 


00. 1 T. 1 1 

Zo:i /:iJ 


rv 'TO 1 rv rvc 
U. /o+U.Uj 


rv rv/i 1 rv rvo 
U.94+U.U0 




1 A A I f\ ^ ^ 
i.44+U.iD 


1 1 1 O'T 

i.io+U.U/ 




00 

u.zz 


ac 
U.Jj 


T/l QOOO 

U. /4J9UZ 


B 


Ti 1 no 1 niiA. 
J i iUZ+ /ZZo 


1 1 .rvo.rv/i 1 
i i.UZ.U4.i 


'TO.O^./IO 

/Z.ZD.4Z 


rv /io_i_rv rv/i 
U.4Z+U.U4 


rv ^c_i_rv rvc 
U.Dj+U.Uj 


00_i_0 1 1 

U.>'Z+U.i i 


1 ^0_i_0 1 

i.oZ+U.iZ 


1 CI _i_o 0^ 
i.Ji+U.Uo 




in 
U.J / 


U.OO 


1 /I ^0000 
i .4DUUUU 


B 


Ti 1 /IT 1 /irvrvrv 
J 1 14/+41KJU 


1 1 . /I T. rvi rv 

i i:4/:Ui.U 


/irv.rvrv.Q/^ 

4U:UU:Jd 


rv crv 1 rv rv/i 
U.J9+U.U4 


rv 00 1 rv rv/^ 
U.oZ+U.Uo 


01 1 1 1 
U.91+U.1 1 


1 01 1 1 
i.Ji+U.iZ 


1 1 /I 1 oc 
i.i4+U.Uj 




CO 

U.JO 


O/I 

U.o4 


1 000000 
i.UooUUU 


B 


TI 1 CO rvo/io 
J i i jZ— Uo4Z 


1 1 .CO. 1^1 
ii.jZ.iO.i 


rvo./io./irv 
—Uo.4Z.4U 


rv TCj^rv rvc 
U. /j+U.Uj 


1 1 Qj^rv rvrv 

i.i j+u.uy 


1 oc j^O 1 
i.Zj+U.Zi 




1/11 j^o 00 
i.4i+U.Uo 


^'Tj^O f\1 

U.D/+U.UJ 


C/^ 
U.JD 


T/1 

U. /4 


Q TOOOO 

Z.J /UUUU 


B 




1 1 . CQ . 1 rv 

ii:jJ:io.U 


4^: JU: 14 


1 1 1 rv rv/^ 
i.lD+U.Uo 


1 CQ 1 rv rv/^ 
l.jJ+U.Uo 


1 /^O 1 1 
l.oU+U.lZ 


1 '71 1 1 1 

1. /J+U.i 1 


1 00 1 oc 
I.oU+U.Uj 




/I /I 

U.44 


1 T 
U.i / 


'J QOO 1 

U.JJJ9oi 


B 


TI Q 1 rv 1 Qoo 1 
J i JiU+jZZi 


1 Q . 1 rv.o^ 
i J. iU.Zo.Z 


Qo.oi .crv 
JZ.Zi. jU 


1 Qrvj^rv rv^ 
i.JU+U.Uo 


1 c^j^rv rvrv 
i.JD+U.Uy 


O'Tj^O 1 

Z.Z/+U. io 


O^j^O 1 ^ 

Z.^D+U. io 


Q coj^o 00 
J.jo+U.Uo 




Q OQ 

J.UJ 


1 /I c 

i.4j 


00^000 

u.yyouuu 


T3 \7 

B, V 


TI QQO rvcrvrv 


1 'J . Q . rvrv 
iJljZlUy.D 


rvc. rvrv. i a 

-Uj:U9:i4 


rv O/I 1 rv rv/i 
U.o4+U.U4 


1 rvrv 1 rv rvrv 
1.UU+U.U9 




1 /^/^ 1 QO 

i.DD+U.JZ 


1 10 1 1 
i.JU+U.iU 


^0 1 00 
U.Do+U.UJ 


U.oJ 


/IT 

U.4/ 


1 coooo 
Z.ijUUUU 


B 


TI QQ'? 1 OC^ 


1 . Af\ 1 
i J.J /.4U.i 


1 O.C/C./l^ 
— iZ. JD.4D 


TTj^rv rv/c 
Z. / /+U.UD 


Q /irvj^rv 1 rv 
J.4U+U.iU 


Q 'TOj^O oq 

J. /y+u.zj 


/I A j^O O/I 

4.40+U.Z4 


A QO j^O 1 

4.JZ+U.iU 


^ o^j^o 00 
D.Uo+U.Zy 


c 00 
j.UU 


O/I 

Z.o4 


CQOOOO 

u.jjyuuu 


T3 \7 
B, V 


T 1 /I C /I 

J1457-J543 


1 /I . cn . ^ A n 

i4:57:z4.7 


-35:43:19 


rv c 1 1 rv rv/i 

0.51+0.04 


rv 'T/C 1 rv rvT 

0.76+0.07 


/I 1 1 

0.84+0.16 


1.59+0.17 


1 1 1 00 
i.io+U.Oo 


00 1 oc 

0.90+O.Oj 


00 1 O/I 

0.00+0.04 


00 1 oc 

0.9j+0.0j 


1 /1 /1 000 
i.4Z4UUU 


B 


J1506+4237 


15:06:52.1 


42:37:26 


0.41+0.04 


0.60+0.07 




1.00+0.18 


0.58+0.06 




0.41 


0.41 


0.587000 


B 


J1516+0014 


15:16:37.0 


00:14:06 


1.00+0.05 


1.12+0.09 


1.34+0.18 


2.29+0.20 


1.71+0.08 




0.96 


1.59 


0.052489 


B 


J1549+0237 


15:49:28.1 


02:37:34 


1.14+0.04 


1.29+0.08 


1.50+0.17 


2.22+0.16 


1.78+0.07 




0.92 


1.11 


0.414414 


B, V 


J1553+1255 


15:53:06.0 


12:55:08 


0.38+0.05 


0.60+0.07 




1.57+0.23 


0.96+0.07 




0.41 


0.74 


1.290000 


B 


J1604+5714 


16:04:45.8 


57:14:35 


0.26+0.04 


0.46+0.06 




1.02+0.07 


0.91+0.03 




0.49 


0.33 


0.720000 


B 


J1635+3807 


16:35:16.3 


38:07:34 


2.66+0.05 


3.51+0.08 


3.86+0.13 


4.57+0.15 


3.96+0.07 




2.40 


3.22 


1.813570 


B 



Table 4. GPS source candidates in the ERCSC — Continued 



Name 


RA 


Dec 


^143 GHz 


^100 GHz 


^ 70 GHz 




^ 44 GHz 


^ 30 GHz 


^ 20 GHz 


^8.6 GHz 


^ 4.8 GHz 


Redshift 


Notes 








[Jvl 


rJvi 


[JyJ 






[Jvl 


rJvi 

vy J i 


[Jvl 


rJvi 

\y J i 








T 1644-77 16 


16:44:32.4 


-77:16:05 


56+0 05 


77+0 08 


0.87+0. 


16 






90+0 08 

yj . y yj — i— yj »yj\j 


40+0 02 

yj . 1 v_/ _!_ yj . yj ^ 


0.24+0.01 


25+0 01 

yj , — !_ yj , yj i 


042700 

yj . yji^ / yjyj 


B 


J 1703-62 13 


17:03:39.4 


-62:13:08 


0.96+0.04 


1.33+0.07 


1.48+0. 


16 


2, 


.67+0.14 


2.24+0.08 


1.05+0.05 


1.03+0.05 


1.04+0.05 




B 


J1727+4530 


17:27:33.6 


45:30:14 


0.75+0.03 


0.94+0.05 


1.10+0. 


11 


1, 


.19+0.14 


1.38+0.05 




1.36 


0.94 


0.717000 


B 


J1927+7358 


19:27:51.6 


73:58:37 


2.42+0.04 


3.28+0.05 


3.82+0. 


10 


4, 


.94+0.11 


5.01+0.05 




3.70 


3.63 


0.302100 


B 


J2005-1821 


20:05:15.6 


-18:21:04 


0.46+0.05 


0.62+0.08 


0.98+0. 


17 






1.18+0.09 


0.60+0.03 


0.42+0.02 


0.40+0.02 


0.868000 




J2006+6424 


20-06-17 8 

m\J\J m ± / m\J 


64:24:07 


36+0 07 


49+0 08 


0.78+0. 


13 


1 


.14+0.11 


1 13+0 06 




0.96 


0.72 


1 574000 

1 . ^ / ^yjyjyj 


B 


J2009+7229 


20-09-58 1 


72-29-06 


49+0 04 


86+0 06 


1.24+0. 


10 


0, 


.90+0.10 


68+0 04 

yj »yj\j — !— yj . v_/ i 




0.79 


0.91 




B 


J2035-6845 


20:35:25.0 


-68:45:04 


0.51+0.06 


0.65+0.08 


0.90+0. 


14 


1, 


.35+0.17 


1.02+0.07 


0.47+0.02 


0.37+0.02 


0.41+0.02 


1.084000 


B 


J2 100-2932 


21:00:57.1 


-29:32:02 


0.51+0.04 


0.71+0.07 


1.02+0. 


15 


1, 


.17+0.21 


0.62+0.10 


0.39+0.02 


0.48+0.02 


0.56+0.03 


1.492000 


B 




21:26:36.5 


-46:07:41 


0.44+0.06 


0.69+0.08 






1, 


.24+0.17 


1.04+0.07 


0.55+0.03 


0.64+0.03 


0.69+0.03 


1.670000 


B 


J2139+1424 


21:39:08.6 


14:24:14 


0.39+0.06 


0.73+0.08 


1.30+0. 


17 


2, 


.48+0.13 


2.64+0.06 




2.27 


1.07 


2.427000 


B 


J2147-7536 


21:47:06.7 


-75:36:40 


2.12+0.04 


2.48+0.06 


3.09+0. 


14 


3, 


.37+0.12 


2.93+0.06 


0.49+0.03 


0.86+0.04 


1.01+0.05 


1.139000 


B 


J2225+2119 


22:25:36.0 


21:19:23 


1.10+0.05 


1.19+0.08 


1.28+0. 


12 


2, 


.16+0.15 


1.87+0.06 




1.40 


1.02 


1.959000 


B 


J2230-3940 


22:30:45.6 


-39:40:19 


0.45+0.06 


0.72+0.07 


1.06+0. 


14 






0.90+0.06 




0.53 


0.56 


0.318049 


B 


J2235-4836 


22:35:01.7 


-48:36:32 


0.62+0.04 


0.93+0.05 


1.07+0. 


12 


1, 


.19+0.14 


1.32+0.06 


1.99+0.09 


2.07+0.10 


1.21+0.06 


0.510000 


B 


J2239-5701 


22:39:16.3 


-57:01:55 


0.80+0.04 


1.02+0.06 


1.33+0. 


13 


2, 


.03+0.13 


1.99+0.06 


0.93+0.05 


1.01+0.05 


0.85+0.04 




B 


J2253+1609 


22:53:59.8 


16:09:07 


27.94+0.16 


28.83+0.14 


27.91+0.29 


22.75+0.20 


16.84+0.10 




10.38 


14.47 


0.859000 


B 



The 143 to 30 GHz flux densities of these sources are from the ERCSC (Planck Collaboration 2011c). The 20 GHz flux density is from the AT20G catalogue (Murphy et al. 2010). If a source 
is in the AT20G catalogue, its 8.6 GHz and 4.8 GHz flux densities are from the AT20G catalogue. If not, the 8.6 flux density (8.4 GHz to be accurate) and 4.8 GHz flux densities are obtained from 
the CRATES catalogue (Healey et al. 2007). Note that the CRATES catalogue does not provide flux uncertainties. The tabulated redshift values are from the NASA/IPAC Extragalactic Database 
(NED). 

In the "Notes" column, "B" indicate sources in the Roma-BZCAT blazar catalogue (Massaro et al. 2008) and CGRaBS blazar catalogue (Healey et al. 2008), "V" means these sources are found 
to be variable in Torniainen et al. (2005). 

Source J2253+1609 is the well known quasar 3C454.3 that only shows GPS-like spectrum during strong flares (Rachen et al. 2010). See also discussion in § 5.2. 
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Fig. 7. SED plots of J0308+0405 (upper) and J2005-1821 
(lower). Planck data are in red, WIMP-Tyr data (if available) 
are in blue, archival data obtained from NED are in green. 



5.1.3. Broad-peaked radio sources 

As noted above, many sources showing GPS-like features are 
radio sources with dominant compact jet emission, usually iden- 
tified with blazars. At VLBI resolution, they generally appear 
as one-sided, compact, jet-like objects. The most likely explana- 
tion for the very flat, sometimes inverted, radio spectra of these 
sources is the overlay of synchrotron- self absorbed emission in 
a continuous jet (IMarscher 1977; ]V[arscher & Gear 1985). We 
may note at this point that the application of this mechanism 
to a standard Blandford-Konigl plasma jet (Blandford & Konigl 
1979) yields an optically thick spectral index a ^0.3-0.5, de- 
pending on the electron spectral index (]V[arscher & Gear 1985). 
The turnover to a completely optically thin spectral index a < 0, 
which is expected to happen somewhere in the GHz-THz regime 
from jet- size considerations, would then naturally produce a 
"GPS type" spectrum. We may therefore conclude that GPS-type 
blazars agree better with the simplest model of compact, contin- 
uous jets than do "typical" blazars which show a low frequency 
spectral index a (Planck Collaboration 201 le). 

IMost of these objects are known to be strongly variable, on 
time scales down to one day (usually on week-month scale in the 
millimeter-regime, however). Variability and one-sided VLBI 
morphology give strong arguments in favor of beaming, i.e., the 
emission is strongly Doppler boosted, making the measured flux 
of the source at a given frequency depend on the Doppler fac- 



tor D as D^~^ . Therefore, small changes in the Doppler boost- 
ing, as expected in helix-shaped or precessing jets, can lead to 
strong variability. The flux density changes would be expected to 
be achromatic, i.e., leaving the spectral shape unchanged. This 
can be distinguished from variability due to distortions in the 
jet, which is expected to be emphasised in a specific frequency 
range; such sources are discussed in § 5.2. 

Figure 8 shows two examples of apparently achromatic vari- 
ability: J1800-h7828 shows a rather typical GPS blazar spectrum, 
with a ^0.3 below the peak frequency at about 10 GHz, and 
steepening to a -0. 17 for y < 100 GHz, with a further break to 
0^ ^ -0.5 at still higher frequencies; The more dome-like spec- 
trum of J0423-0120, with a ^ 1 at lower frequencies, points to a 
more uniform component producing the radiation. Supporting, 
ground-based, radio data from Eff'elsberg, IMetsahovi and the 
VLA further suggest significant, achromatic variability for this 
source, as would be expected from changes of the Doppler boost- 
ing in a helical or precessing jet. The high frequency spectrum 
has the same characteristics as J1800-h7828. The temporarily 
very flat index a ^ -0.2 after the peak may suggest an unusu- 
ally flat electron spectrum, as discussed in Planck Collaboration 
(201 le). 



5.2. Flat and multi-component spectra 

In our examination of radio sources in the ERCSC we found 
a large number with extended, flat, power law spectra. For in- 
stance, the flux density of J2203-hl725, as determined by Planck 
and VLA observations, varies only between 0.99 and 1.18 Jy 
between frequencies of 4.9 and 100 GHz, and decreases only 
to 0.906 Jy at 143 GHz. IMany of these flat spectrum sources 
are the blazars discussed in Planck Collaboration (201 le). We 
also found several sources with prominently zig-zag or "bumpy" 
spectra. In some cases the variation in flux density from one 
Planck band to the next was several times the associated errors. 
This could result from the superposition of emission from sev- 
eral components, as is probably the case in the examples shown 
in Figure 9. However, it is important to recall that PlancJ^s multi- 
frequency observations were not exactly simultaneous. A source 
could sweep through the Planck beams at one frequency days 
before doing so at a neighboring frequency (depending on posi- 
tion, it takes 7-10 days for a source to drift entirely across the 
focal plane). Hence fast variability on time-scales of days can 
contribute to irregularities in its Planck spectrum. We expect this 
efl'ect to be smaller for the HFI frequencies, given the relatively 
closer packing in the focal plane. We discuss below the possible 
physical conditions for true multi-component sources, followed 
by examples of spectral artifacts caused by variability. 



5.2.1 . Flat spectrum sources with a variable component 

Spectra resulting from superposed radiation of at least two com- 
ponents have become the standard model to explain flaring 
blazars; this is known as a "shock-in-jet" model (Marscher & 
Gear 1985). One component, the jet, produces a flat (a~0) 
spectrum up to some break frequency, mostly between ~10 and 
-100 GHz, above which the spectrum steepens to a typical in- 
dex a — 0.7. The second component, associated with an evolv- 
ing shock, produces a self-absorbed synchrotron spectrum with 
a turnover frequency ~ 100 GHz. Below this turnover frequency, 
the spectrum is strongly inverted. As the shock evolves, the 
turnover frequency and the flux density of the shock component 
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Fig. 8. VLA, IMetsahovi, IRAIM, Effelsberg and Planck measurements of two Icnown GPS blazars: Jl 800+7828 (left) and J0423- 
0120 (right). While the former shows only small signs of variability, the latter is clearly variable while preserving its spectral shape, 
pointing to helical or pressing jet (see text). Asterisks identify epochs within 10 days from Planck observations at any of the LFI 
channels. For both sources, the Planck ERCSC data are a superposition of two scans separated by 3 and 6 months, respectively. 



changes, and in some phases a detectable contribution of this 
component may not be present at all. 

An example of this class is J2253+1609, better Icnown as 
3C454.3. During an outburst in 2005/2006, the source showed a 
strong spectral pealc around 300 GHz, evolving to lower frequen- 
cies over several months, while observations made in 2004 did 
not show this feature (see Villata et al. 2007, Raiteri et al. 2008, 
Villata et al. 2009, Rachen et al. 2010 and references therein). 
Incidentally, this source showed a similar outburst when it was 
first observed with Planck (Planck Collaboration 201 le), which 
is why it is included in Table 4. Two other examples where 
Planck observations play a crucial role are shown in Figure 9. 
The sources shown in Figure 5, all identified with known blazars, 
present a very similar spectral shape. We therefore expect that 
their high-frequency-peaked spectrum is a temporary feature. 

We emphasise that Planck, which is expected to perform at 
least four full sky surveys, is a powerful tool to distinguish vari- 
able from permanent spectral features. By comparing catalogues 
derived from individual surveys made by Planck, we will be able 
to resolve the questions raised here in the full frequency regime 
accessible to Planck, even without reference to external data. 
Results from this research will be presented in future Planck pa- 
pers. 

5.2.2. Artifacts in the ERCSC spectra of variable sources 

In this subsection, we discuss some potential problems in in- 
terpreting ERCSC spectra for variable sources. As noted in § 1, 
44 GHz observations of a given source take place at two times 
separated by 7-10 days, for each scan. If the source is strongly 
variable on a time scale of order days, the 44 GHz flux density 
contained in the ERCSC could be an awkward average. An ex- 
ample is the well-known, rapidly variable sources J0722-h7120 
(Ostorero et al. 2006), shown in Figure 10. Supporting observa- 
tions, particularly the IRAIM data at 86 and 143 GHz, show it to 
have varied on short time scales during the Planck survey. The 
apparently anomalous 44 GHz Planck measurement is an aver- 
age of the flux density in the source's high and low luminosity 
state. We show this case as a warning: spectral anomalies seen 
only at 44 GHz may be due to this efl'ect. 



In sources with still faster variability, the efl'ect of non- 
simultanity of the Planck observations can produce even more 
bizarre efl'ects. The ERCSC spectrum of the source Jl 159+2914, 
also known as the IDV blazar TON599, shows a strong zig- 
zag shape, dropping by a factor of two between 44 and 70 
GHz, followed by another small bump (see Figure 11). This 
source is known to show very fast variability, and the compar- 
ison with Efl'elsberg data, and IRAIM 86 GHz and 143 GHz data 
taken around the time of the Planck scans, suggests that it had 
a strong flare at a peak frequency of -50 GHz in the first days 
of June 2010. This flare must have started after 23 IMay (when 
the Efl'elsberg observations were made), and probably declined 
again on 7 June when the Planck 100 GHz point was taken. The 
ERCSC spectrum of the source is a superposition of a quite low 
state in December 2009, and the high, flaring state in June 2010, 
except for 30 and 70 GHz measurements, which were made by 
Planck for the second time shortly after the last day of data used 
for ERCSC, and are therefore not included in the average. 

These examples make clear that a fair amount of care has to 
be taken when intepreting Planck ERCSC spectra for variable 
sources. 



5.3. Sources showing evidence of a spectral upturn at high 
frequency 

In many radio sources (Arp 220 as an example), synchrotron 
emission is dominated at high frequencies by re-emission from 
warm dust. As already noted, essentially no extragalactic Planck 
sources show this pattern at frequencies below 143 GHz. Of all 
the sources at | | > 30° that have Planck measurements at 30 to 
217 GHz, only eleven show evidence of a significant increase of 
flux density even from 143 GHz to 217 GHz; one of these is the 
nearby and well- studied star-forming galaxy NGC253 (Figure 
12). Of the remaining 10 sources, 7 are HII regions located in the 
]V[agellanic clouds, like J0047-7310 and J0048-7306; two are 
Galactic sources; and the remaining one is 1M82, a well-known 
starburst galaxy. The lack of upturn- spectrum sources in the 
ERCSC suggests that most of the sources detected by Planck are 
flat spectrum sources like blazars, with high enough synchrotron 
luminosity to swamp dust reemission even up to 217 GHz. 
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Fig. 9. Two examples of clearly multi-component spectra. Planck measurements are in filled circles; AT20G data are in open 
diamonds; PACO data, with their dense sampling in frequency space, are shown in open squares. Asterisks identify epochs within 
10 days from Planck observations. The multi- valued spectra in J2035-6846 are a clear sign that this source is variable on time scales 
less than two weeks (see text). 
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Fig. 10. Based on ground-based and Planck data, the apparently anomalous 44 GHz point in the spectrum of the strongly variable 
source J0722-h7120 is shown to be due to the averaging of flux densities described in § 5.2.2. Asterisks identify epochs within 10 
days from Planck observations at any of the LFI channels. 



We also note here the potential danger of apparently upturn- 
ing spectra being produced by source confusion. The low proba- 
bility of source confusion discussed in § 4. 1 is valid only for the 
low frequency catalogues; at 545 and 857 GHz, in particular, the 
ERCSC catalogues are much richer, and these frequencies are in 
the Ray leigh- Jeans part of the spectra for most IR sources. The 
apparently upturning spectra may result from confusion of a ra- 
dio source disappearing in the background with an IR source ap- 
pearing within the large beam of Planck Careful checks with in- 
frared catalogues are required to verify any case of a potentially 
upturning spectrum. Hence we restrict our attention in this paper 
to frequencies of 30-353 GHz, and generally to 30-217 GHz. 

5.4. ERCSC LFI sources with no plausible match in existing 
radio catalogues 

Finally, we looked for sources with such extreme spectra that 
they have no counterparts in existing low-frequency, large area, 
radio surveys. This serves the dual purpose of looking for un- 
usual or extreme sources and validating the ERCSC. We looked 
for matches in the following catalogues: the WIMAP 7 year 



catalogue (Gold et al. 2011), the NEWPS catalogue (IMassardi 
et al. 2009), the 20 GHz AT20G survey in the southern hemi- 
sphere (IMurphy et al. 2010), the CRATES catalogue in the north- 
ern hemisphere (Healey et al. 2007), and the GB6 catalogue 
(Condon et al. 1994). To reduce the number of spurious iden- 
tifications of the bright Planck sources, a flux density cut of 
0.3 Jy was applied to both the AT20G and GB6 catalogues. In 
the northern hemisphere, where the high frequency coverage 
is incomplete, we also looked at the 1.4 GHz NVSS catalogue 
(Condon et al. 1998) with a 0.5 Jy flux density cut. We cross- 
correlated these catalogues with the ERCSC at 30 GHz, 44 GHz 
and 70 GHz using a search radius equal to 0.5 FWHIM at each 
Planck frequency channel. We have shown in § 4 that this search 
radius is unlikely to yield spurious matches. We also experi- 
mented with difl'erent flux cuts, but found that our adopted cuts 
are reasonable compromises between matching most sources 
and avoiding spurious matches. Significantly lower flux cuts re- 
sult in spurious matches because the chance of random associa- 
tion given the large LFI beams goes up. The unmatched sources 
listed in Table 5 for example often have NVSS sources within the 
Planck LFI beams, but the 1 .4 GHz flux densities of these NVSS 
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Fig. 11. The zig-zag shape spectrum of the strongly variable source Jl 159+2914. Asterisks identify epochs within 10 days from 
Planck observations at any of the LFI channels. 
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Fig. 12. Planck sources with an upturn in the spectrum. Left: J0047-7310 in the SMC. Center: Planck spectrum of a familiar 
star-forming system, NGC253. Right: The starburst galaxy M82. 



sources are on the order of a few mJy; thus they are clearly 
unlikely to be associated with the ERCSC sources, whose flux 
denisities are on the order of 1 Jy or above. 

With these automatic matching steps, we found that among 
high Galactic latitude (\b\> 5°) sources, 12 ERCSC sources at 
30 GHz, 13 at 44 GHz and 26 at 70 GHz do not match with 
any known source. These represent respectively 2 %, 5 % and 
8 % of the high Galactic latitude sources in each of these cata- 
logues. However, there is a caveat in such automatic matching 
procedure. While the above matching radius is reasonable as- 
suming that Planck dominates the positional uncertainty, this is 
not necessarily the case, for example, when matching against 
the WMAP sources. The comparatively smaller beam sizes at 44 
and 70 GHz of Planck suggest that this procedure could miss 
genuine WMAP matches. 

Since sources that exist in only one Planck channel are more 
likely spurious, we first looked for unmatched sources that ex- 
ist in more than one Planck band, and found a number of such 
cases among the 44 and 70 GHz "unmatched" sources. We then 
looked in NED for each of these unmatched sources to see if 
they have counterparts in other radio catalogues that are not 
among those listed above. In particular, extended sources may 
appear in older, lower resolution surveys such as PMN (Griffith 
& Wright 1993), but not in newer, higher resolution surveys such 
as AT20G. Several such cases were found. We further associated 
a number of unmatched sources with Galactic objects such as 
PNe and SNR. After these checks, the number of residual, ap- 
parently unmatched sources is 5 at 30 GHz, 6 at 44 GHz, and 
15 at 70 GHz. These are listed in Table 5. Judging from the 



Galactic latitude, 2 of the 30 GHz, 3 of the 44 GHz and 4 of the 
70 GHz sources are likely associated with the Galaxy while one 
of the 70 GHz sources is likely in the LMC; these are flagged 
in Table 5. The remaining sources are either exciting "new" 
sources or else spurious. To distinguish between these possibil- 
ities, we looked at the postage stamp images for all the sources 
that remain unidentified (see Figure 13). In all cases, the postage 
stamps are 5 FWHM on a side. A number of the sources at 44 
and 70 GHz appear to be of low significance and hence could 
be spurious. Some are flagged in the ERCSC as possibly con- 
taminated by CMB signals. Others are at low Galactic latitude. 
One apparent exception at 70 GHz is G158. 33-20. 53 (but see 
comments below). At 30 GHz, some sources are of high signifi- 
cance but appear extended, which, given the ~30' Planck beam 
size, means they are likely to be either Milky Way objects or as- 
sociated with very nearby galaxies (or could be CMB artifacts). 
An exception is G2 19. 17-08.90 (5*30 ghz = 1 Jy) which appears 
to be real and point- like. From NED, we find that there is a 
^60 ;um = 68 Jy IRAS source (IRAS06282-0935) near this posi- 
tion that seems to have been identified as IR cirrus (Strauss et al. 
1992). The most convincing 70 GHz source, G158.33-20.53, 
also has a counterpart, IRAS03259+3105, in NED, which again 
is identified as cirrus by Strauss et al. (1992). 

We thus conclude that we do not identify any genuinely new 
population of sources among the ERCSC catalogues at the LFI 
frequencies. While confirming a match is generally much eas- 
ier than claiming a source is "unmatched", we can say that af- 
ter the automated procedure described above, > 90 % of the LFI 
sources had a reasonable match in an external radio catalogue. 
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Fig. 13. ERCSC sources with no plausible identifications in low frequency radio catalogues. 



After cross-matching between the three LFI bands and search- 
ing in NED for plausible identifications, especially important in 
the case of extended radio sources, we are left with a total of 26 
potentially "unidentified" LFI sources, almost all of which can 
probably be explained by artifacts or extended Galactic struc- 
tures. 



6. Conclusions 

We summarise in this section the primary conclusions to be 
drawn from Planck's study of extragalactic radio sources. We 
stick closely to the observational results, and provide comments 



on the fit between these observational results and current theo- 
ries of the physics of extragalactic radio sources. 

6.1. Overall properties 

Planck has demonstrated that the high frequency counts (at least 
for frequencies < 143 GHz) of extragalactic sources are domi- 
nated at the bright end by synchrotron emitters, not dusty galax- 
ies. This finding is in agreement with conclusions reached by the 
South Pole Telescope (SPT) and Atacama Cosmology Telescope 
(ACT) teams (Vieira et al. 2010; Marriage et al. 2010, respec- 
tively) from the properties of sources at much lower flux den- 
sities. An inference from this result is that the cores of extra- 
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Table 5. ERCSC sources at | /? | > 5° with no obvious match in 
external radio catalogues 



ERCSC Name 


Notes* 


PLCKERC030 G219. 17-08.90 


G 


PLCKERC030 G252.09+46.04 




PLCKERC030 G27L21-09.56 


G 


PLCKERC030 G328.86-8L59 




PLCKERC030 G354.78+23.68 


G 


PLCKERC044 G014.94+29.84 




PLCKERC044 G075.52-08.69 


G 


PLCKERC044 G105.35+09.85 


G 


PLCKERC044 G105.43-07.07 


G 


PLCKERC044 Gl 8 1.40+56.02 




PLCKERC044 G359.05-23.36 




PLCKERC070 GOOO.86-27.44 


G? 


PLCKERC070 G030.76+05.22 


G 


PLCKERC070 G03 1.58+05.37 


G? 


PLCKERC070 G053.99-10.45 


G 


PLCKERC070 G095.46+45.89 




PLCKERC070 G120.02+09.88 


G 


PLCKERC070 G158.33-20.53 


G? 


PLCKERC070 G212. 11-19.20 




PLCKERC070 G235.56-51.54 




PLCKERC070 G245.36-20.14 




PLCKERC070 G250.22+24.04 




PLCKERC070 G279.96-33.42 


L? 


PLCKERC070 G292.76-40.83 




PLCKERC070 G339.56+09.30 


G 


PLCKERC070 G353.92+ 15.84 


G? 



* Here G stands for likely Galactic source. L stands for likely LMC 
source. 



galactic sources, and not extended structure, dominate the high 
frequency emission, as suggested by earlier work at frequencies 
below most of the Planck bands (see recent discussions in Lin 
et al. 2009 and Murphy et al. 2010). The conclusion that the 
core dominates the high frequency emission is supported by the 
close agreement between Planck and VLA flux densities (§4), 
implying that the emission is found within the VLA beam of 
order arcseconds in size. 

The emerging dominance of emission from a flat spectrum 
core implies that extrapolation of flux densities and/or counts 
of sources to frequencies > 30 GHz cannot be reliably made 
from low frequency catalogues, where the emission is domi- 
nated by lobes. This has been recognised for some time (see Lin 
et al. (2009) and Murphy et al. (2010)); Planck strongly con- 
firms this conclusion. If we look at the SEDs of the sources, 
including all components, the spectra grow flatter as frequency 
increases and the steep- spectrum lobes fade away. On the other 
hand, Planck observations allow us to follow the SEDs of a sta- 
tistically significant sample of sources to much higher frequen- 
cies than can generally be employed in ground-based observa- 
tions. As noted in § 4 (and described in more detail in Planck 
Collaboration 201 Id), Planck provides clear evidence of a spec- 
tral steepening in the radio and millimeter wave emission from 
extragalactic cores; the steepening sets in at frequencies above 
-44-70 GHz. The observed spectral steepening in turn means 
that radio sources contribute less foreground noise to increas- 
ingly sensitive searches for small angular scale anisotropics in 
the CMB than earlier models had suggested (e.g., de Zotti et al. 
2005). This finding, too, is consonant with results from SPT and 
ACT. 



6.2. Properties of individual sources 

Planck allows us for the first time to investigate sub-mm spec- 
tra beyond 200 GHz, which is the limit of most ground-based 
monitoring programs. In this regime, the spectra of radio sources 
are usually expected to be optically thin synchrotron emission, 
and our results largely confirm that this part of the spectrum 
is represented by a single power law. The vast majority of ex- 
tragalactic sources in the ERCSC lists at 30-100 GHz are flat 
spectrum radio sources of the sort discussed in detail in Planck 
Collaboration (201 le), with a scattering of bright steep spectrum 
sources strong enough to register in the lowest frequency bands 
of Planck. A small fraction (< 10%; Table 4) show peaked or 
convex spectra, but even these sources are mostly identified as 
known blazars. As blazars are known to be variable, care has to 
be taken in the interpretation of their spectra. For most sources, 
ERCSC flux densities represent the average over two scans, and 
even during one scan, spectral artifacts may occur due to the fact 
that it requires 7-10 days for the Planck focal plane to cross a 
point source. 

In contrast, most known examples of compact, newborn ra- 
dio galaxies (CSOs), which were originally thought to be the 
dominant class of sources with Gigahertz-peaked spectra, are 
mostly too faint to be detected by Planck. A new population of 
bright, very compact, high-peaked CSOs has not been found. 
Likewise, very few of the extragalactic radio sources found by 
Planck show evidence of the sharp spectral upturn expected from 
dust reemission at high frequencies. Only NGC253 shows a clear 
upturn at a frequency < 143 GHz; sources that show dust reemis- 
sion dominating at 217 GHz are in many cases Galactic sources, 
in the Magellanic clouds, or nearby known star-forming galax- 
ies. 

Although we have not investigated in detail all the sources 
with no obvious matches in lower frequency radio catalogues 
(§ 5.4), let alone every extragalactic source in the ERCSC, we 
find no convincing evidence for the emergence of a new and un- 
expected population of sources. The 26 ERCSC sources with no 
match in radio catalogues appear to be a heterogeneous mixture 
of conventional radio sources, many of them Galactic. 

6.3. Future observations and analysis 

Two sets of future observations will help clarify the status of 
some of the sources listed in § 5.1 and § 5.2. Careful monitor- 
ing and/or VLBI observations of sources with peaked spectra 
would determine whether they are highly compact radio galaxies 
(CSOs), compact jet objects with stable GPS type spectra, or are 
instead merely flaring sources with temporarily convex spectra. 
In fact, our results seem to suggest most extreme spectral fea- 
tures seen in Planck sources may be associated with flares. With 
at least four full sky surveys, Planck is a powerful tool to distin- 
guish these possibilities, by deriving catalogues from individual, 
six-month, sky surveys and comparing flux densities of all bright 
point sources at six-month intervals. Further ground-based ob- 
servations, now underway at frequencies below and overlapping 
the Planck frequency bands, will support this efl'ort. 

Likewise, for the non-variable extreme radio sources like 
CSOs, the addition of more surveys to the Planck maps will al- 
low us to extract deeper catalogues, with the potential to find 
some more of these usually faint objects. 

In the case of Galactic sources, the role of CO emission lines 
is important. The situation for extragalactic sources is more com- 
plicated, because the CO lines redshift in and out of the Planck 
bands. We have underway a study to look for the influence of 
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CO emission in the ERCSC spectra. We expect the effect to be 
small, since most of the ERCSC sources are extremely bright, 
non-thermal emitters, with strong continua. There is no obvious 
evidence in the SEDs we have examined for the presence of CO 
emission. In this regard, we again warn readers to be careful in 
interpreting anomalous 44 GHz observations; a spectral bump at 
44 GHz is not necessarily evidence for redshifted CO. 

We used ground-based observation data for both validation 
purposes, specifically cross-calibration between ground-based 
instruments and Planck and in the study of the spectral and vari- 
ability properties of extreme radio sources. Most of these obser- 
vation programs, like the Planck mission itself, are still ongoing. 
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